[a,b,x,y]=solve('a+b=1','a*x+b*y=2/3','a*x^2+b*y^2=2/5','a*x^3+b*y^3=2/7','a','b','x','y');